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Abstract 

A  procedure  based  on  dynamic  fracture  mechanics  is  proposed  for  assessing  the 
effectiveness  of  tear  straps  in  a  rupturing  airplane  fuselage  weakened  by  a  row  of 
multiple  site  damage  (MSD).  A  large  deformation,  elastic-plastic  finite  element  model  of 
the  rupturing  fuselage  with  an  unsymmetrical  crack  flap  is  used  to  demonstrate  the 
existence  of  a  mixed  mode  I  and  II  crack  tip  deformation  and  a  large  axial  stress 
preceding  the  propagating  crack.  These  numerical  results  were  used  to  evaluate  the 
dynamic  crack  curving  and  crack  arrest  criteria  and  hence  to  assess  the  possibility  of 
crack  curving  as  the  crack  approaches  a  tear  strap  without  the  presence  of  MSD. 


Introduction 

The  role  of  a  tear  strap  in  an  airplane  fuselage  is  to  arrest  an  axial  crack  which  is 

propagating  either  subcritically  under  fatigue  loading  or  dynamically  after  reaching 

criticality.  Such  crack  arrest  can  occur  by  the  lower  stresses  due  to  the  reinforcing 

effect  of  the  tear  strap  or  by  crack  curving  due  to  the  complex  crack  tip  stress  field 

generated  by  the  crack  flap  and  the  tear  strap.  The  latter  provides  controlled  damage  and 
depressurization  of  the  fuselage  Although  the  tear  strap  is  the  last  defense  of  an  axially 
rupturing  fuselage,  little  is  known  about  its  effectiveness  in  arresting  an  axial  crack  and 
much  of  its  design  is  based  on  empirical  rules  derived  from  sub-  and  full-scale  testings 
of  pressurized  fuselages.  Literature  is  abundant  with  analytical  and  experimental  papers 
dealing  with  axial  fatigue  crack  extension  and  possible  arrest  in  idealized  and  actual 
fuselages  but  few  consider  crack  curving  and  flapping  as  an  arrest  criterion.  The  only 
analytical  paper  dealing  with  tear  straps  and  crack  curving  and  flapping  as  an  arres! 
criterion  appears  to  be  that  of  Kosai  and  Kobayashi  (1991). 

In  the  above  mentioned  paper,  the  axial  crack  was  assumed  to  open  symmetrically,  i.e. 
mode  I  crack  tip  deformation,  in  a  fuselage  reinforced  by  longitudinal  stringers,  frames 
and  tear  straps.  If  failure  were  to  occur  along  the  multiple  site  damages  (MSD's)  in  a  lap 
joint,  evidences  [Sampath  and  Broek  (1991);  NTSB  (1989))  show  that  fracture  will  occur 
only  in  the  upper  skin  with  the  lower  skin  still  attached  to  the  fasteners  The  axial 
stretching  due  to  the  resultant  one-sided  crack  flap  would  then  generate  axial  tensile 
stresses  along  the  cracked  upper  skin  and  impose  a  mode  II  state  of  crack  tip  deformation 
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The  importance  of  mode  II  deformation  and  the  attendant  mode  II  stress  intensity  factor. 
K||,  was  also  observed  recently  in  a  fatigue  crack  growth  study  of  a  miniaturized  and 
idealized  model  of  a  stringer  reinforced  fuselage  by  Fyfe  and  Sethi  (1991). 

The  purpose  of  this  study  is  to  refine  the  finite  element  modelling  of  axial  crack 
propagation  and  curving  in  a  pressurized  fuselage  [Kosai  and  Kobayashi  (1991)]  by 
accounting  for  this  more  realistic  mode  II  crack  tip  deformation. 


Finite  Element  Modelling 

Due  to  the  unsymmetric  crack  tip  deformation  of  an  axial  crack  along  a  lap  joint,  one  half 
of  the  airplane  fuselage,  as  shown  in  Figure  1,  must  be  modelled  by  finite  elements. 
Implicit  in  this  modelling  is  an  assumption,  which  was  made  to  reduce  the  computing 
time,  that  an  identical  axial  crack  was  also  propagating  in  the  other  symmetric  half  of 
the  fuselage.  The  geometric  parameters  associated  with  the  fuselage  and  the  material 
properties  of  the  2024-T3  clad  aluminum  used  in  the  skin  and  the  tear  strap  are  shown  in 
Table  1. 

One  of  the  unknown  quantity  is  the  actual  pressure  distribution  on  as  well  as  in  the 
vicinity  of  the  crack  flap  as  it  propagates.  As  shown  in  Figure  2,  in  one  analysis,  the 
crack  flap  was  assumed  to  be  loaded  with  the  cabin  pressure  through  the  layer  of  cabin 
insulation  which  acted  as  a  weak  bladder  and  which  maintains  the  cabin  pressure  during 
the  short  subsecond  duration  of  rapid  axial  crack  propagation.  In  another  analysis,  the 
effect  of  a  leakage  in  cabin  air  was  studied  by  prescribing  a  linearly  decreasing  pressure 
on  the  crack  flap  while  maintaining  the  cabin  pressure  on  the  uncracked  portion  of  the 
fuselage. 


Figure  3  shows  the  finite  element  model  (FEM)  of  the  rupturing  fuselage.  The  longitudinal 
stringers  were  modelled  by  beam  elements,  the  frames  were  replaced  with  beam  elements 
and  the  straps  were  represented  by  increased  skin  thickness.  The  limit  cabin  pressure 
differential  of  51.7  kPa  (7.5  psi)  was  prescribed.  A  quasi-static,  elastic-plastic  finite 
element  analysis  was  used  to  model  the  rapid  crack  propagation  in  this  ductile  material. 
This  quai-static  analysis  was  justified  since  the  measured  crack  velocities  in  the  thicker 
and  more  brittle  7075-T6  and  7178-T6  aluminum  plates  was  less  than  ten  percent  of  the 
dilatational  wave  velocity  [Kobayashi  and  Engstrom  (1967);  Arakawa,  Drinnon,  Kosai  and 
Kobayashi  (1991)]. 

Also  the  limit  hoop  stress  in  the  uncracked  skin  of  the  fuselage  is  a  low  1.06  MPa 
(15.4  ksi)  and  thus  the  computed  crack  tip  plastic  zone  size  is  less  than  one  of  the  coarse 
mesh  element  when  the  crack  length  less  than  a  ■  178  mm  (7  in.).  Thus,  the  prior  history 
of  plastic  yielding  for  a  crack  length  less  than  178  mm  (7  in.)  was  ignored  and  an 
incremental  static,  elastic-plastic  finite  element  analysis  was  initiated  at  a  *  178  mm 
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(7  in.)  by  advancing  the  crack  tip  one  nodal  distance,  i.e.  25.4  mm  (1  in.),  and  by 
superposing  the  prior  elastic-plastic  state  to  that  after  each  crack  advance. 

The  crack  curving  criterion,  which  is  triggered  as  the  crack  approaches  the  tear  strap  is 
based  on  the  near  crack  tip  state  of  stress  and  thus  considerable  refinement  of  the  finite 
element  mesh  was  necessary  to  meet  the  required  numerical  accuracy.  Refined  mesh,  on 
the  other  hand,  increased  exponentially  the  computation  time  of  this  incremental  elastic- 
plastic  analysis.  In  order  to  reduce  the  computer  time,  the  crack  tip  region,  where  crack 
kinking  was  likely  to  occur,  was  first  identified  through  a  coarse  mesh  finite  element 
analysis.  This  region  was  then  reanalyzed  with  a  refined  mesh  using  a  rezoning  procedure 
as  shown  in  Figure  3.  This  simulation  of  crack  propagation  was  continued  by 
incrementally  advancing  the  crack  tip  using  the  coarse  mesh  model  and  by  repeating  the 
analysis  with  the  refined  mesh  using  the  restart  file  generated  from  the  coarse  mesh 
model.  The  crack  kinking  angle  was  then  computed  and  the  kinked  crack  path  was 

generated  incrementally  by  using  the  refined  mesh  in  the  crack  tip  region  using  the 

restart  file  generated  from  previous  analysis. 

Crack  Curving 

As  mentioned  previously,  the  function  of  a  tear  strap  is  not  so  much  as  to  arrest  a 
propagating  axial  crack  by  reducing  the  circumferential  stress  in  the  crack  path,  but  to 
deflect  the  crack  in  the  circumferential  direction.  The  large  opening  due  to  crack  flapping 
would  then  reduce  the  crack  driving  force  through  controlled  depressurization  of  the  cabin 
and  thus  arrest  the  crack.  Such  crack  deflection  is  accomplished  by  the  presence  of  mode 
II  stress  intensity  factor,  K||,  which  is  generated  by  the  flap  of  the  upper  skin.  The  large 

crack  flap  also  generates  a  large  axial  stress  ahead  of  the  propagating  crack  tip 

[Kobayashi,  Emery,  Love  and  Chao  (1988);  Kobayashi,  Emery,  Love  and  Chao  (1988)],  and 
together  with  the  reduced  circumferential  stress  due  to  the  presence  of  the  tear  strap 
augments  the  propensity  for  crack  curving  under  mixed  mode  crack  tip  loading. 

In  a  previous  analysis  (Kosai  and  Kobayashi  (1991)],  an  elastic  crack  kinking  criter'on  for 
a  stationary  mode  I  crack  [9,10]  with  plasticity  correction  was  used  to  assess  the 
effectiveness  of  the  tear  strap  for  a  symmetric  crack  flap.  The  unsymmetrir  srack  flap 
considered  in  this  study,  however,  will  generate  a  combined  mode  I  anc  II  crack  tip 
deformation  field  [Fyfe  and  Sethi  (1991)]  and  thus  crack  kinking  is  inherent  in  such  crack 
tip  deformation. 

If  the  apparent  mode  I  stress  intensity  factor,  K|,  is  elevated  by  tHe  presence  of  a  row  of 
MSD’s  along  the  axis  of  the  crack,  then  the  crack  will  continue  to  propagate  in  the  axial 
direction  despite  the  presence  of  an  inevitable  crack  tip  bulging  and  the  resultant  K||.in  a 
pressurized  fuselage  [NTSB  (1989)].  The  physical  evidence  of  such  self  similar  crack 
extension  along  a  row  of  MSD’s  in  the  presence  of  Kn  thus  requires  a  new  crack  curving 
criterion. 
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Using  the  maximum  stress  criterion  with  the  crack  tip  stress  field,  the  angle  of  crack 
kinking,  9c  in  the  presence  of  Kn  was  derived  by  Erdogan  and  Sih  (1963).  This  maximizing 
condition  yields 


K|  sin  9c  +  K,i  (3  cos  9c  -1)  =  0  (1) 

where  K|  and  Kn  are  the  modes  I  and  II  stress  intensity  factors  respectively. 

This  crack  curving  criterion  predicts  a  positive  9c  for  a  negative  Kn  since  K|  is  always 
positive  as  shown  in  Figure  4.  Likewise,  a  negative  9c  is  predicted  for  a  positive  Kn.  In 
the  absence  of  Kn,  however,  Equation  (1)  predicts  a  self  similar  crack  propagation  or  9c  = 
0  and  fails  to  explain  the  physically  observed  crack  instability  where  crack  curving  takes 
place  in  a  K|  field.  The  static  elastic  crack  curving  criterion  of  Streit  and  Finnie  (1980) 
or  Ramuiu  and  Kobayashi  (1983)  predicts  such  crack  curving  and  was  found  to  agree  well 
with  available  experimental  data.  In  the  following,  a  brief  account  of  a  mode  II  extension 
of  this  crack  curving  theory  is  given. 

The  mixed  mode,  elastic  crack  tip  stress  field  in  terms  of  the  crack  tip  polar  coordinate 
system  is  given  by 


1  9  9  3  9  a0x 

arr  *  -j=z  cos  g  [K|(1  +  sin2  ^  ^  Kn  sin  9-2  Kn  tan  -  ]  +  g  (1  +  cos  29) 

1  9  9  3  C[ox 

006  COS  2  [KIC0S2  2  -  2  K|1  Sin  0]  +  2  H'005  2°)  (2) 

1  9  CTox 

Tre  -  —  cos  g  [Ki  sin  9  +  Kn  (3  cos  9  -  1)]  -  g  sin  20 


where  a0x  is  the  second  order  term  and  is  commonly  referred  to  as  the  remote  stress 
component.  The  same  maximum  circumferential  stress  criterion  is  then  used  to  derive  a 
mixed  mode  crack  extension  criterion  by  assuming  that  fracture  will  occur  when  the 
maximum  circumferential  stress  is  equal  to  the  equivalent  circumferential  stress  in 
mode  I  fracture.  Experimental  evidences  compiled  by  Ramuiu  and  Kobayashi  (1983) 
heavily  favors  the  incorporation  of  the  second  order  term,  i.e.  Equation  (3)  in  predicting 
the  crack  kinking  angle  under  mixed  mode  fracture.  The  mixed  mode  fracture  criterion  of 
Equation  (1)  then  becomes 


Kic 


„  *  9c 

K|  cos3  Y 


„  9 c  9c 
3Kn  cos2  y  sin  Y  ♦ 


a0x  (1  -  cos  29) 


(3) 
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Equation  (3)  incorporates  the  second  order  effect  into  Equation  (1)  but  otherwise  follows 
the  general  trend  of  Equation  (1).  If  Equation  (1)  predicts  a  positive  8c  for  a  positive  K| 
and  a  negative  K||,  Equation  (3)  shows  that  for  the  same  K|  and  Ku  and  positive  a0x,  0c  =  0 
and  the  crack  will  propagate  straight  ahead.  If  K|  >  Ku,  Equation  (3)  also  shows  that  for 
the  same  Ku  and  c0x.  crack  kinking  is  only  possible  when  K|  <  Kic. 

The  crack  kinking  angle  can  be  obtained  by  maximizing  the  circumferential  stress  and 
results  in  the  following  transcendental  equation: 


Kji 

K| 


-sin 


ec 


(3  cos  6c- 1 ) 


[2  cos 


16  V  2  it 


A*  COS  8C] 


(4) 


A  / - °OX 

A=Vrc  (5) 

The  term  A,  is  related  to  the  critical  distance,  rc,  from  the  crack  tip  and  is  proportional 
to  the  nonsingular  stress,  oox  Ramulu  and  Kobayashi  (1983)  has  shown  that  rc  is  a 
material  dependent  parameter  which  must  be  determined  experimentally.  Thus,  unlike  the 
crack  curving  criterion  represented  by  Equation  (3),  the  crack  kinking  angle  computed  by 
Equation  (4)  incorporates  the  second  order  term  and  represents  the  mixed  mode  extension 
of  crack  kinking  criterion  of  Streit  and  Finnie  (1980)  and  Ramulu  and  Kobayashi  (1983). 
As  shown  in  Figure  5,  the  crack  kinking  angle,  0c.  increases  with  increasing  aox  and 
negative  a0x  tends  to  stabilize  the  crack  path. 

The  crack  extension  criterion  represented  by  Equation  (3)  and  the  crack  kinking  criterion 
represented  by  Equation  (4)  does  not  account  for  the  elevation  in  stress  intensity  factors, 
K|  and  possibly  Ku,  due  to  the  presence  of  MSD.  The  effect  of  MSD  obviously  is  to  promote 
self-similar  crack  extension,  as  seen  in  the  NTSB  report  (1989),  and  can  be  incorporated 
into  this  analysis  by  artificially  increasing  K|  by  a  magnification  factor  which  represents 
the  interaction  effect  between  the  axial  crack  and  the  small  crack  emanating  from  the 
adjacent  fastener  hole. 

Finally  in  the  absence  of  an  elastic-plastic  crack  curving  criterion,  the  above  elastic 
crack  curving  criterion  was  used  to  estimate  the  crack  curving  potential  of  a  tear  strap. 


Results 

Unlike  the  fuselage  considered  in  (Kosai  and  Kobayashi  (1991)],  the  crack  flap  is 
constrained  from  opening  by  the  adjacent  stringer  which  is  213  mm  (8.4  in.)  away  from 
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the  extending  axial  crack.  The  smaller  flap  opening  in  turn  resulted  in  a  small  plastic 
yield  zone  at  the  crack  tip  as  alluded  to  in  previous  sections 


As  mentioned  earlier,  an  elastic  analysis  was  first  conducted  to  study  the  extent  of  crack 
tip  plastic  yielding  in  order  to  determine  the  crack  length  at  which  incremental  elastic- 
plastic  analysis  would  be  used  with  each  increment  of  crack  extension.  This  incremental 
analysis  was  initiated  when  the  crack  tip  yield  region  reached  the  size  of  the  crack  tip 
element  or  at  a  crack  length  of  178  mm  (7  in.).  This  elastic  plastic  analysis  was 
conducted  until  the  crack  length  reached  a  =  508  mm  (20  in.).  The  resultant  K|,  Ku  for  the 
limit  pressure  differential  of  51.7  kPa  (7.5  psi)  computed  by  the  coarse  and  fine  grid 
finite  element  analyses  are  shown  in  Figures  6  and  Figure  7,  respectively.  Also  shown  in 
these  figures  are  K|  and  Ku  without  the  tear  strap  (TS),  but  with  the  stringer  and  frame 
This  FEM  analysis  models  either  a  fuselage  without  a  tear  strap,  or  a  fuselage  with  a 
disbonded  tear  strap. 

Figure  6  shows  that  the  first  crack  kinking  occurred  at  25.4  mm  (1  inch)  before  the  tear 
strap.  The  subsequent  kinked  crack  analysis  used  a  refined  grid  for  each  incremental 
crack  extension  while  superposing  successive  numerical  results  on  previous  analyses. 
Figure  7  shows  that  the  kinked  crack  propagated  in  a  straight  line  after  which  it  kinked 
again  very  close  to  the  tear  strap.  While  the  predicted  kinking  angle  without  the  non¬ 
singular  stress,  i.e.  using  Equation  (1),  is  about  0c  *  30  degrees,  the  predicted  kinking 
angle  with  the  non-singular  stress,  i.e  using  Equation  (4),  is  about  0C  *  40  degrees.  Figure 
6  also  shows  little  differences  in  K|  and  Ku  with,  and  without,  the  tear  straps. 

Figure  9  shows  the  crack  path  obtained  through  this  analysis.  The  results  appear  to  agree 
qualitatively  with  those  observed  experimentally  [Maclin  (1991)].  Also  shown  in  Figure  8 
is  the  finite  element  idealization  of  the  skin-stringer-frame  construction. 

Figure  8  shows  the  resultant  K|,  K||  when  a  linearly  varying  cabin  pressure,  as  shown  in 
Figure  2,  was  prescribed  on  the  crack  flap.  A  comparison  of  Figures  6  and  9  shows 
negligible  differences  between  the  stress  intensity  factors  associated  with  the  different 
flap  pressure.  Thus  details  in  the  pressure  distribution  in  the  vicinity  of  the  axial  crack 
may  be  unimportant  as  long  as  the  cabin  pressure  is  maintained  through  the  bulk  of  the 
fuselage  during  crack  propagation. 

The  flap  deformations  due  to  the  two  pressure  boundary  conditions  are  shown  in  Figure  2. 
The  fully  pressurized  flap  deformed  9  mm  (0.35  in.)  more  than  the  partially  pressurized 
flap  with  little  differences  in  the  stress  intensity  factors  as  mentioned  above. 


Conclusions 

The  utility  of  large  deformation,  elastic-plastic  finite  element  analysis  with  rezoning  for 
predicting  crack  curving  in  a  pressurized  fuselage  with  tear  straps  was  demonstrated.  An 
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elastic  crack  propagation  and  kinking  criteria  in  the  presence  of  mixed  mode  fracture 
were  presented. 


Discussion 

The  large  K||  and  the  large  axial  stress  ahead  of  the  propagating  crack  were  developed 
through  the  large  unsymmetric  crack  flap.  For  an  airplane  in  flight,  the  high  velocity  air 
flow  along  the  fuselage  will  tend  to  close  the  flap  thus  counteracting  the  effect  of  the 
opening  flap  due  to  cabin  pressure.  The  net  result  is  that  the  axial  crack  will  continue  to 
propagate  without  kinking  and  hence  without  arresting.  Perhaps  the  large  opening  in  the 
Aloha  Airlines  B737-200  [NSTB  (1980)]  could  have  been  due  in  part  to  this  air  flow  in 
flight  which  suppressed  the  development  of  crack  flap  necessary  for  crack  curving  and 
hence  arrest. 
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Table  1 .  Structural  and  material  property  data  for  aircraft  fuselage  analysis 


Tear  Strao 


Figure  1 


Figure  2 


Axial  rupture  of  an  idealized  fuselage. 
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lower  Skin  — 

Section  A- A 

Constant  prmura  on  flap 
(  a-457.2  ) 


Stringer 


Section  A-A 
Line  arty  varying  praaaura  on  flip 
(  1-457.2  ) 


Estimated  pressure  on  crack  flap. 
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A 


Tear  Strap  (Shdl  Element) 
Frame  (Beam  Element) 


Coarse  Grid  FEM  Model 


Figure  3.  FEM  meshing  procedure 


Mode  II 


Mode  II 


Figure  4.  Crack  extension  angle  under  mixed  mode  loading. 
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Frame 


Hall  Crack  Length,  a  (mm) 

Figure  7.  Stress  intensity  factor  variations  with  crack  extension  (fine  grid  analysis), 
p  =  51.7  kPa  (7.5  psi). 


Figure  8.  Crack  path  near  tear  strap  and  finite  element  idealization  of  skin-stringer- 
frame  construction. 
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